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Reply by A.N.Ivanov and P. Kienle: Nowadays it is well established experimentally that neutrinos i^a with 
Icpton flavours a — e,^ and r are superpositions = J2j ^ajl^j) rnassive neutrino mass-eigenstates {vj) with 
masses ruj, where U*^ are elements of the 3x3 unitary mixing matrix U, defined by mixing angles 0ij [H- The wave 
functions {va) and {vj) are orthogonal and used for the description of neutrino oscillations < — ^ with frequencies 
ujij = Amfj/2£', where E is the neutrino energy and Amf^- = mf — In K-shell electron capture {EC) decays 

of the H~like heavy ions m — > d + i^e, where m and d are mother and daughter ions in their ground states 
one deals with an emission of electron neutrinos |z^e) = Ylj Thus, the i^C-decay rates of the H~like heavy 

ions are defined by the decay channels m — >■ + Vj, where the final states are described by the orthogonal wave 
functions {vidi\djVj) = for i 7^ j. The states of the daughter ions dj differ in 3~momenta qj and energies Ed{qj). 
The massive neutrinos Vj are produced with 3~momenta kj ~ —qj and energies Ej{kj), caused by conservation of 
energy and momentum in the decay channels m ^ dj + Vj . Since massive neutrinos Vj are not detected they appear 
. . . in the asymptotic states with 3-momenta kj, energies and energy differences uJij = Amfj/2Mm, where is 

" the mass of the mother ion m In the GSI experiments [2, Si the i?C-decay channels m — >■ dj + Vj are measured 
■ by detecting the daughter ions dj. If the daughter ions would be detected in the asymptotic states with 3-momenta 
[ qj and energies Ediqj), the probability per unit time of the i?C-decay m — > d + i^e is equal to 

; P{m^ dv,){t) = ^ \U,,\^P{m ^ d, v,){t) = |[/,,f d,iy,m\\ (1) 

• where A{m — > djh'j)(t) is the amplitude of the decay channel m dj + Vj. However, this is not the case in 
\ the GSI experiments, where the time differential detection of the daughter ions with a time resolution leads to 
indistinguishability of daughter ions in the decay channels m — >■ d j + Vj . As a result the daughter ions dj are measured 
in the asymptotic state d with a 3-momentum q and an energy Ed{q) such that q — qj and Ed{q) — Ed(qj) 01- This 
does not violate the orthogonality of the wave functions in the final state {vid\dvj) = for i ^ j. The energy and 
^ _ momentum uncertainties 5Ed and \5qd\, respectively, induced by the time differential detection of the daughter ions, 
^ • provide the overlap of the wave functions of the daughter ions \i 5Ed ^ | ujij \ and \5qd\ S> \qi~qj\ — \ki — kj\, where 
P5 \ present also the differences of the recoil energies of the daughter ions. The time differential detection of the daughter 
ions dj in the asymptotic state with the 3-momentum q and an energy Ed{q) results in a smearing of momenta and 
^sj ■ energies in the decay channels m — dj + Vj around q + kj ~ and Ed{q) + Ej{kj) ~ M„j. This is the origin of the 
^ ' non-vanishing interference terms in the probability per unit time of the i?C-decay m ^ d + Ue Q 
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P{m dve){t) = \Uej\^P{m -> dvj){t) + 2 ^ — Re[C/;,f7ej.4*(m dv,){t)A{m dvj){t)l (2) 
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. where in comparison with Eq.(l) the second sum is caused by the interference terms. Since uncertainties SEd and 
I \Sqd\ are rather small 0] and \kj \ = \qj\ ~ \q\ ~ Qec and Qec ^ f^j, where Qec is the Q-value of the i?C-decay 
■ m ^ d + Ve, one can set the neutrino masses zero everywhere except for energy differences ujij and mixing angles 9ij 
in the interference terms 0], and neutrino 3-momenta equal kj ~ k. As a result the -EC-decay rate is given by Q 

" 1 /" d^a d^k / \ 

X '--^'^ - 2M^y ^ (2.)32i^, (2.)32i;. - + 2^Re[l/* t/e,]cos(..,0)- (3) 

Thus, the asymptotic orthogonality of the final state wave functions does not influence the observation of the time 
modulation of the -EC-decays, observed in the GSI experiments with a time resolution <^ Tij much shorter than 
the modulation periods Tij = 2'!r /cuij This is unlike the assertion by Flambaum Q. For Td 3> Tij the daughter 
ions dj become distinguishable with 3-momenta_g^ and energies Ed{qj) the time modulation vanishes (see Eq.(l)). 
For the theoretical description of the GSI data y, accounting for the procedure for the detection of the daughter 
ions, one can use time-dependent perturbation theory and wave packets for the wave functions of the daughter ions, 
related to the density matrix description of unisolated quantum systems Q . For technical details we refer to 0] . 
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